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SECTION  I 
INTRODUCTION 

The  use  of  transmission-line  theory  plays  an  important  role  in 
the  analysis  of  system  response  to  a nuclear  electromagnetic  pulse  (EMP) . 
Various  transmission-line  models  are  typically  employed  for  EMP  internal 
Interaction  studies  so  as  to  calculate  how  EMP  energy  is  distributed 
within  the  interior  regions  of  the  system. 

Inputs  or  source  terms  for  such  transmission-line  models  are 
series  voltage  generators  and  shunt  current  generators  distributed  along 
various  portions  of  the  transmission  lines.  Such  sources  arise  from 
EMP  energy  penetrating  through  slits,  holes,  or  other  apertures  in  the 
outer  shielding  surface  of  the  system,  or  diffusing  through  skin  panels 
and  then  coupling  to  the  transmission  line. 

The  response  of  interest  is  usually  a voltage  or  current  at  a 
particular  point  on  the  line,  usually  at  an  impedance  element  which 
represents  the  load  impedance  of  a circuit  which  might  be  damaged  or 
upset  by  the  EMP  transient.  In  many  cases,  both  time  haraonic  and 
transient  responses  are  considered  to  permit  the  prediction  of  EMP 
effects  on  the  system. 

One  of  the  most  straightforward  approaches  for  solving  the 
internal  interaction  problem  involves  modeling  a complicated  trans- 
mission-line network  by  one  or  more  single-wire  transmission  lines  which 
are  interconnected  in  scries  or  parallel,  but  with  no  closed  loops 
being  formed.  Such  single-line  models  are  described  in  references  (1) 
and  (2).  In  addition  to  these,  references  (3)  through  (9)  discuss 

1.  Plonsey,  R. , and  R.E.  Collin,  Principles  and  Applications  of 
Electromagnetic  Fields,  McGraw  Hill,  1961. 

2.  Tesche,  F.M. , et  al..  Evaluation  of  Present  Internal  EMP  Interaction 
Technology;  Description  of  Needed  Improvements,  AFWL-TR-75-288,  Air 
Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  NM,  Oct.  1973. 

3.  Tesche,  F.M. , and  T.K.  Liu,  Selected  Topics  in  Transmission-Line 
Theory  for  EMP  Internal  Interaction  Problems,  AFWL-TR-77-73,  Air 
Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  NM,  Dec.  1976. 
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various  improvements  which  can  be  made  to  the  single-line  model  by  more 
accurately  accounting  for  the  nonidealized  geometry  found  in  the  actual 
interaction  problem. 

Another  approach  for  treating  the  internal  interaction  problem 
is  to  use  multiconductor  transmission-line  theory.  While  this  Increases 
the  complexity  of  the  calculations  involved  for  a particular  problem, 
the  accuracy  to  which  the  computational  model  approximates  the  actual 
problem  is  much  greater. 

The  anlaysis  of  multiconductor  transmission  lines  is  well 
documented  in  the  literature.  One  popular  approach  is  to  employ  a 
lumped  parameter  network  (LPN)  model  for  the  multi conductor  transmission 
line  and  then  perform  a circuit  analysis  using  a large  network  analysis 
code.  Such  a procedure  is  discussed  in  reference  (10) . A major 
difficulty  with  such  an  approach  lies  in  the  fact  that  if  a complicated, 
electrically  long  cable  harness  is  to  be  analyzed,  substantial  time  must 
be  spent  in  defining  the  network  to  be  analyzed,  choosing  the  element  values, 
and  setting  up  the  problem.  Moreover,  computation  time  can  be  substantial. 

4.  Tesche,  F.M. , and  T.K.  Liu,  An  Electric  Model  for  a Cable  Clamp 
on  a Single-Wire  Transmission  Line,  AFWL-TR- 76-325 , Air  Force 
Weapons  Laboratory,  Kirtland  Air  Force  Base,  NM,  Mar.  1977. 

Coen,  S.,  T.K.  Liu,  and  F.M.  Tesche,  Calculation  of  the  Equivalent 
Capacitance  of  a Rib  near  a Single-Wire  Transmission  Line,  AFWL-TR- 
77-60,  Air  Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  NM,  June  1977. 

6.  Lee,  K.S.H.,  and  F.C.  Yang,  A Wire  Passing  by  a Circular  Aperture 
in  an  Infinite  Ground  Plane,  AFWL-TR-77-52,  Air  Force  Weapons 
Laboratory,  Kirtland  Air  Force  Base,  NM,  Feb.  1977. 

7.  Lam,  J.,  Equivalent  Lumped  Parameters  for  a Bend  in  a Two-Wire 
Transmission  Line,  AFWL-TR-77-5,  Air  Force  Weapons  Laboratory, 

Kirtland  Air  Force  Base,  NM,  Jan.  1977. 

8.  Lam,  J.,  Propagation  Characteristics  of  a Periodically  Loaded 
Transmission  Line.  AFWL-TR-76-324 , Air  Force  Weapons  Laboratory, 

Kirtland  Air  Force  Base,  NM,  Jan.  1977. 

9.  Liu,  T.K.,  Electromagnetic  Coupling  between  Multiconductor  Trans- 
mission Lines  in  a Homogeneous  Medium,  AFWL-TR- 7 6- 33 3,  Air  Force 
Weapons  Laboratory,  Kirtland  Air  Force  Base,  NM,  May  1977. 


10.  EMP  Handbook  for  Missiles  and  Aircraft  in  Flight.  AFWL-TR-73-68,  Air 
Force  Weapons  Laboratory,  Kirtland  Air  Force  Base,  NM.  Also,  EMP 
Interaction  Notes,  Vol.  1-1. 
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An  alternate  approach  for  treating  the  multiconductor  transmission 
line  is  to  solve  the  set  of  transmission-line  equations  directly.  The 
solution  is  accomplished  by  using  the  modal  concept,  as  discussed  in 
references  (11)  through  (18) . This  method  involves  determining  the 
various  quasi-TEM  modes  and  related  propagation  constants  which  can 
exist  on  the  multiconductor  transmission  line.  This  is  achieved  by 
determining  the  eigenvalues  and  eigenvectors  of  a propagation  matrix 
which  contains  the  fundamental  line  parameters.  The  total  solution 
for  voltage  or  current  on  the  multiconductor  line  is  then  descrlbable 
as  a linear  combination  of  all  possible  modes  on  the  line. 

Previous  application  of  multiconductor  transmission-line  theory 
to  the  Internal  interaction  problem  has  Involved  simple  branching  trans- 
mission lines  or,  at  times,  a single  section  of  line.  To  solve  the 
transmission-line  network  problem  with  branching,  all  branches  except 
one  are  represented  by  a Thevenin  or  Norton  equivalent  circuit  and  the 


11.  Carson,  J.R.,  "The  Rigorous  and  Approximate  Theories  of  Electrical 
Transmission  along  Wires,"  Bell  System  Technical  Journal,  No.  1, 

Jan.  1928,  pp.  11-23. 

12.  Kuznetsov,  P.I.,  and  R.L.  Stratonovich,  The  Propagation  of  Electro- 
magnetic Waves  in  Multiconductor  Transmission  Lines,  Pergamon  Press, 
Oxford,  England,  1964. 

13.  Kajfez,  D.,  "Multiconductor  Transmission  Lines,"  AFWL  EMP  Interaction 
No tes , Note  151,  June  1972. 

14.  Paul,  C.R.,  "On  Uniform  Multimode  Transmission  Lines,"  IEEE  Trans. 
M.T.T. , Vol.  MTT-21,  No.  8,  Aug.  1973,  pp.  556-558. 

15.  Paul,  C.R.,  "Efficient  Numerical  Computation  of  the  Frequency 
Response  of  Cables  Illuminated  by  an  Electromagnetic  Field,"  IEEE 
Trans.  M.T.T. , Vol.  MTT-52,  No.  4,  Apr.  1974,  pp.  454-457. 

16.  Paul,  C.R.,  "Useful  Matrix  Chain  Parameter  Identities  for  the 
Analysis  of  Multiconductor  Transmission  Lines,"  IEEE  Trans.  M.T.T., 
Sept.  1975,  pp.  756-760. 

17.  Frankel,  S.,  Cable  and  Multiconductor  Transmission  Analysis,  Harry 
Diamond  Laboratory,  HDL-TR-091-1,  June  1974. 

18.  Lenahan,  T.A.,  "The  Theory  of  Uniform  Cables,  Parts  I & II," 

Bell  System  Technical  Journal,  April  1977,  pp.  597-625. 
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resulting  single-section  transmission  line  is  solved  to  obtain 
voltages  and  currents.  This  process  is  repeated  until  all  branches  in 
the  network  are  all  studied.  Not  only  is  this  method  tedious  and 
time  consuming,  but  it  has  the  weakness  of  not  being  able  to  treat 
transmission-line  networks  with  complicated  branching  or  connections 
involving  closed  loops. 

In  many  practical  cases,  however,  the  wiring  configuration  is 
not  simple.  In  general,  there  may  be  a variety  of  closed  loops  formed 
by  transmission  lines  in  an  actual  system,  and  Induced  currents  on  this 
transmission-line  network  can  behave  quite  differently  from  currents 
induced  on  a simple  branching  transmission-line  network. 

The  method  adopted  in  this  report  of  analyzing  a general  lossless 
multiconductor  transmission-line  network  is  based  on  the  scattering 
matrix  technique  used  widely  in  distributed  network  analysis  (ref.  19). 

Thus,  at  a junction  where  a number  of  transmission-line  tubes  meet,  the 
reflected  currents  (or  voltages)  at  the  nodes  are  related  by  a scattering 
matrix  to  the  incident  currents  (or  voltages)  at  all  the  nodes  at  the 
same  junction.  These  quantities  at  the  two  ends  of  a branch  are  also 
related  by  a propagation  matrix.  With  the  knowledge  of  excitation 
sources  and  termination  conditions,  these  two  relations  are  sufficient 
to  yield  the  current  (and  voltage)  at  all  nodes  by  means  of  matrix 
inversion.  This  procedure  is  similar  to  the  solution  of  lumped  circuit 
problems  using  mesh  analysis  (ref.  20),  and  has  been  discussed  recently 
in  the  literature  (refs.  21  through  23). 

19.  Collin,  R.E.,  Foundations  for  Microwave  Engineering,  McGraw  Hill, 

New  York,  1966. 

20.  Desoer,  C.A.,  and  E.S.  Kuh,  Basic  Circuit  Theory.  McGraw  Hill, 

New  York,  1969. 

21.  Baum,  C.E.,  "The  Role  of  Scattering  Theory  in  Electromagnetic  Inter- 
ference Problems,"  in  Electromagnetic  Scattering.  P.L.E.  Uslenghi  (ed.), 
to  be  published. 

22.  Baum,  C.E.,  "Coupling  into  Coaxial  Cables  from  Currents  and  Charges 
on  the  Exterior,"  presented  at  the  1976  URSI  Meeting.  Amherst,  Mass. 

23.  Tesche,  F.M. , "A  General  Multiconductor  Transmission-Line  Model." 
presented  at  the  1976  URSI  Meeting.  Amherst,  Mass. 
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This  report  is  intended  to  give  a concise  summary  of  the 
solution  procedui  s for  lossless  multiconductor  transmission-line 
networks.  The  detail  description  of  the  basic  theory  is  discussed  in 
a separate  report  (ref.  24).  Reference  (24)  also  treats  the  more  general 
transmission-line  networks  that  may  include  losses. 

Various  definitions  and  concepts  are  Introduced  in  Section  II. 
The  overall  procedure  of  solving  the  transmission-line  network  problem 
is  also  summarized.  In  Section  III,  an  example  of  how  to  apply  this 
method,  as  well  as  some  numerical  results,  are  given. 


24.  Baum,  C.E.,  T.K.  Liu  a^d  F.M.  Tesche,  "On  the  General  Analysis  of 
Multiconductor  Transmission-Line  Networks,"  AFWL  EMP  Interaction 
Notes,  to  be  published. 
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SECTION  II 

BASIC  CONCEPT?  AND  APPROACH 


In  this  section,  a few  definitions  and  concepts  that  are 
unique  to  general  transmission-line  networks  are  first  introduced. 

These  quantities  are  useful  for  the  discussions  in  the  later  part  of 
this  section  and  the  rest  of  the  report.  Detailed  descriptions  of 
these  quantities  are  presented  in  reference  (24).  The  fundamental 
relations  needed  for  the  development  of  the  computer  code  are  then 
summarized. 

1.  Graph  Representation 

In  Figure  la,  a simple  multiconductor  transmission-line  network 
is  shown.  In  analogy  to  the  lumped  network  analysis,  a graph  is  intro- 
duced in  Figure  lb  to  represent  the  entire  network.  The  basic  elements 
of  the  graph  are  tubes  and  junctions.  A tube  is  a section  of  .:ulti- 
conductor  transmission  line  with  uniform  or  gradually  changing  configura- 
tions. In  this  report,  we  assume  that  there  is  one  common  reference 
conductor  for  the  whole  network. 

A junction  is  where  one  or  more  tubes  meet.  We  denote  a 

junction  by  , j=l,2,...,Nj  , where  Nj  is  the  total  number  of 

junctions  within  the  entire  network.  A tube  connected  between 

junctions  J,  and  J,  is  denoted  by  T . . Thus,  T , Identifies 

the  same  tube  as  T,  . , assuming  there  is  only  one  tube  between  the 

J 

two  junctions. 

The  network  graph  representation  enables  one  to  visualize  the 
overall  network  connections  in  a simple  form  and  is  particularly  useful 
when  the  network  is  extremely  complicated  and  contains  a number  of 
closed  loops  and  crossed  connections.  Here,  a double  line  represents  a 
tube  which  contains  many  conductors  and  is  a distributed  element, 
and  a circle  represents  a junction  within  which  multiple  connections 
between  conductors  of  the  tubes  are  possible. 

Figure  lb  does  not  completely  describe  the  network  topology. 
Details  of  Intertube  connections  within  a junction  must  be  described. 
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Junction 


Figure  1.  (a)  A simple  mutliconductor  transmission-line 

network,  and 

(b)  its  associated  graph. 
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A junction  subgraph  is  needed  for  each  junction.  Thus,  in  Figure  2, 
junction  subgraphs  for  the  four  junctions  of  Figure  1 are  presented. 

In  the  junction  subgraph  where  more  details  are  shown,  each 
line  represents  a conductor  in  the  bundle  of  multiconductor  transmission 
lines,  and  each  dot  represents  a node  or  connection  between  individual 
wires  or  impedances.  The  identification  with  single  (or  scalar) 
quantities  thus  makes  it  a more  detailed  representation  compared  with 
the  network  graph  of  Figure  lb,  which  is  identified  with  multiple 
(or  vector)  quantities. 

2.  The  Transmission-Line  Propagation  Matrix 

Current  and  voltage  waves  on  a transmission  line  can  be  decomposed 

into  a forward  traveling  component  and  a backward  traveling  component. 

At  a junction,  these  components  are  also  Identified  as  the  incident 

and  reflected  components.  Let  us  consider  a typical  single  tube  (see 

Figure  3)  which  is  spanned  between  two  junctions  with  voltage  and 

current  sources  at  z = 5 • Here  we  denote  the  two  components  by 

(j(re) (^^0,s))  and  (1^^^ (z=0,s))  at  one  end  (the  left  end)  of  the 

tube  and  (z=L,s))  and  (?'■  \z=L,s))  at  the  other  end  (the 

n n 

right  end).  The  superscripts  (re)  and  (i)  denote  the  reflected  and 
incident  components,  respectively,  at  the  junctions.  The  conditions 
z = 0 or  z = L define  the  junctions  of  interest.  The  parameter  s 
is  the  complex  frequency,  and  L is  the  length  of  the  tube. 

At  the  two  ends  of  the  tube,  the  total  current  is  the  sum  of  the 
two  components  at  that  point,  i.e., 

(f  (z,s))  - (I^^^(z,s))  + (1) 

n n n 


where  z = 0 or  z = L . 

For  a tube  of  multiconductor  transmission  line  with  N 

c 

conductors  plus  a reference,  there  are  N propagating  modes  (refs.  3 

c 

and  14).  In  relating  currents  at  different  positions  along  the  tube,  it 
is  necessary  to  (i)  decompose  the  total  current  into  modal  currents. 
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Figure  2.  Junction  subgraph  at  junction  (a)  1, 
(b)  2,  (c)  3,  (d)  4. 
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(11)  relate  the  nodal  currents  at  different  positions  by  the  modal 
propagation  matrix,  and  (111)  convert  the  modal  components  so  related 
back  to  the  total  current. 


Details  of  this  procedure  have  been  described  In  references  (3) 

and  (24).  Here,  we  summarize  the  Important  results  In  the  following. 

For  the  lossless  case,  the  modal  current  (1  (s))  Is  defined  so  that 

"n 

(t  (z,s))  - (^^  )•(!  (s))  (2) 

n n,m  c^ 

where  ) Is  the  similarity  transform  matrix  formed  by  the 

n ,m 

eigenvectors  of  the  matrix  ((C  )*(L'  ))  . The  corresponding  eigen- 

values  of  (<^  ) are  the  elements  of  the  diagonal  matrix  (y  (s))  , 

n.m 

which  Is  given  by 


(Y^  (s))^  - s^(^„  J -(C:  )*(L’  )•(»^7'  ) 

c n,m  n,m  n,m  n,m 

n.m 


(3) 


The  elements  of  (Y  (s))  are  In  the  form  of  s/(mode  velocity).  The 
c ^ 

characteristic  Impedadce  matrix  (Z  ) Is  given  by 

n,m 

(Z  ) = (C'  )"^(^„  J-(Y  (s))-(^„  )"^/s 

c n,m  n,m  c n,m 

n,m  n,m 


(4) 


The  reflected  current  components  at  the  two  ends  of  the  tube  are  related 
to  the  Incident  components  by 


'(!<"> /cr^w.s)) 


.(i) 


,(!<-> a.s))/  Vf„_,(s))  / yr-a.,)), 


■(i), 


+ (y®^(s)) 


(5) 


where  (p  (s))  Is  the  current  propagation  matrix  of  the  tube  and  Is 
n,m 

given  by 
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(6) 


(f  (s))  - ^ 

n,m  n,m  n.m 

(^^®\s))  is  the  source  vector  of  the  tube  referred  to  the  two  ends. 


and  is  given  by 


(■5‘f  (a))- 


n,in 


•[(li®^5,s))  + (Z„  )"^*(^„^®^5,s))] 


n.m 


-1. /?r(s), 


[(l^®^C.s))  - (Z  )"-"*(^^“'(C.s))] 


\S/ 

If  the  tube  is  excited  by  distributed  voltage  source  (z,s)) 


and  distributed  current  source  (I^  (x.s))  , then 


(.y^^ha))  - j 


L -1 


0 \t  Hexp(Y  (s))[L-5])-(^„  )"^ 

'2  ''n.m  ^ c^^^  n.m 


[(T^®^’(5,s))  + (Z  (s))'^‘(V^^’(C.s))] 
" ‘'n.m 


[(I^®^’(C.3))  - (Z^  (s))"^-(^„^®^’(C.s))] 

n C n 

n.m 


For  the  tube  with  conductors  (plus  a reference), 

is  an  N xN  matrix  and  (3^®2(s))  is  a 2N  xl  vector, 
c c n c 
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At  the  two  ends  of  the  tube  the  total  voltage  Is  given  by 


- q(Z^  )-((r''  ^z.s))  - (r^^z.s))) 

n,m 


where  q * +1  at  z «•  0 and  q ••  -1  at  z ■ L . 


3.  The  Junction  Scattering  Matrix 

For  a junction  with  tubes  connected,  the  quantities 

associated  with  tube  T , are  denoted  by  the  subscript  T , , 


1 = 1,2,. ..,N^  . For  example,  (I  (z,s))_  is  the  current  vector 

k 

of  Tube  T,  , . J’  i 

For  distributed  circuits,  the  reflected  and  Incident  quantities 

at  a junction  are  often  related  by  a scattering  matrix  (ref.  19). 

For  junction  J.  the  current  junction  scattering  matrix  (3  (s)) , 

j n,m  j 

is  defined  so  that  ^ 


. ,s)) 


j J 


V ^’V  i.K 

J, 


where  z , , i * 1,2,...N  is  taken  to  be  the  appropriate  coordinate 

j.ki  J. 

of  the  tube  T , at  the  junction  J.  . 

J 9 J -w 

From  Equation  (9),  it  is  clear  that  the  submatrix  (S 

’ J 

is  given  by  j 


(Z  ,3))  - (S  (S)l.(l^^>(2  S)) 

J**^n  n,m  J.  j.K 

J m 


• -I  . 


Thus,  for  n ■ m,  (S  ),  is  the  current  reflection  matrix, 

n,m  J j 

whereas  for  n m (S„  . is  the  transmission  matrix. 

n,nr 

Two  cases  of  particular  Interest  are  treated  here:  (!)  junctions 
contain  tubes  only  which  are  interconnected,  and  (ii)  junctions  where 
a tube  is  terminated  by  external  Impedances. 

(a)  Junctions  with  tubes  only 

At  a junction  where  there  are  only  direct  connections  be- 
tween tubes,  the  Klrchhoff  current  law  and  the  Kirchhoff  voltage  law 
have  to  be  enforced. 

Kirchhoff 's  current  law  states  that  the  sum  of  the  current 
flowing  out  a node  is  zero.  For  the  case  that  the  n^-th  wire  of  tube 

T.  , is  connected  to  the  n_-th  wire  of  tube  T , , and  to  the 

J , * J » *^2 

, etc.,  then 


j.k^ 


j.k. 


(10) 


t,  • 2-1  t ‘Ij.k  * k ' k ’ 


"j.k 


j.k 


j .k^  j ,k^ 


ly  2,  f ■••yN-  • 


Equation  (10)  can  be  put  into  the  matrix  form,  l.e.. 
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ITK— ^ -M-.-.  TC-.  ^ ■ n mn-jr  - ■ 


i. 


i 


(0  0 ...  q,  . ...  0,0...  q 


I I ^ 

^1...  'o  ... 


^'1  ^j.ki 

ci(z.  . ,s)) 


''j  “j 


In  Equation  (11) , all  elejents  in  the  left  matrix  are  zero  unless 

they  correspond  to  the  conductors  which  are  connected  at  the  node. 

For  Kj  connections  within  the  junction  , there  are  Kj 

equations  similar  to  Equation  (10),  and  we  can  define  the  junction 

connection  matrix  (Cj  )j  so  that 

n,m  j 


(C^  ) • . ^'‘"2 

n.m  I : 

“^j 


n,m 


where  (C_  ) is  a K xM  matrix,  and  M is  the  total  number 

n.m  j j j j 

of  conductors  connected  to  junction  , i.e., 


"j.  ■ Y.  M 


where  N 


Is  the  number  of  conductors  of  tube  T. 


Kirchhoff's  voltage  law,  for  the  case  of  simple  connections, 
requires  all  voltages  associated  with  each  conductor  to  be  the  same 
at  the  same  node.  Thus,  in  the  particular  example  above,  we  have 


V 1,  -r  ~ ^^4  If  T’ 

“l  ^’*^1  ^j,k^  "2  ^j,k2 


\ k T 
“2  ^'^2  ^j,k„ 


k.  'T 


The  above  equation  can  be  easily  written  in  matrix  form.  For  a consistent 
set  of  connections,  there  are  - Kj,  equations.  Let  us  denote  the 


corresponding  matrix  as  (C^  )j^  sucA  that 

n,m  j 


n,m  j 


(V  (z  ,s)) 

" J’^1  ^j.k^ 

(V  (z  ,s)) 

“ J’’"2  ^j,k2 


J 


n,m 


In  addition  to  the  above,  two  special  situations  have  to  be 
considered:  (i)  a conductor  is  open-circuited,  and  (ii)  a conductor 

is  short-circuited  to  ground. 
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is  open-circuited 


When  the  n^-th  conductor  of  tube 
the  current  there  is  zero  and  we  have 


0 


(14) 


This  corresponds  to  the  entry  of  the  value  1 at  the  appropriate 

location  in  the  matrix  (C  ) in  a row  where  all  other  elements 

^ J ■ 

are  zero.  In  this  case,  thet§  is  no  corresponding  equation  for  the 
voltage  law. 

When  a conductor  is  short-circuited  to  ground,  the  voltage  at 
the  node  is  zero,  i.e., 


(15) 


This  corresponds  to  a value  of  1 for  the  appropriate  element  in  the 

matrix  (C^  )j  in  a row  where  all  other _ elements  are  zero. 

n,m 

Using  Equations  (1),  (8),  (12)  and  (13),  the  current  scattering 
matrix  is  given  by 


n,m  J 


n,m  j 

(C^  )j  ‘(qj  k^^c  . 

n,m  j n,m  j,k' 


n,m  j 

n,m  j n,m  j,k 


(16) 
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(b)  Junctions  with  termination  impedance 

For  a tube  terminated  in  an  impedance  network  (Z^  (s))  , 

n,m 

the  current  scattering  matrix  (S^  (s))  relates  the  reflected  and 

,m 

incident  current  components  of  the  tube,  i.e, 

- (§  „(s))-(l^^^z,s))  (17) 

n n,m  n 

In  this  case,  (S  (s))  is  identical  to  the  current  reflection  matrix 
n,m 

(I*  (s))  , which  is  given  by  (ref.  3) 

n 

(S  _(s))  » (F  (s))  - - ((Z^  (s))  + (Z^  ))“^- 


((Z_  (s))  - (Z  )) 

T C 

n,m  n,m 


4.  Solution  of  the  Network  Problem — The  BLT  Equation 

After  evaluating  the  current  propagation  matrices  (P  _(s)) 

(s) 

and  the  source  vectors  (s))  for  all  tubes,  and  the  current 

n 

scattering  matrices  (§  (s)),  for  all  junctions,  one  can  combine 

n ,m  J , 

these  quantities  and  obtain  a network  equation,  which  has  been  called 
the  BLT  equation. 

Each  of  these  above  matrices  or  vector  must  be  arranged  to 
fit  into  the  overall  network  matrices  or  vectors.  This  re-arrangement 
can  be  achieved  by  either  a matrix  transformation  (ref.  23)  or  by 
carefully  renumbering  the  individual  matrix  and  vector  elements  to 
fit  into  the  network  matrices  and  vectors.  Another  way  is  to  consider 
the  network  quantities  as  aupermatrices  or  tensors  of  rank  four  (ref.  24) 
For  the  network,  we  have 


and 


(i^“\s))  - (P  (S)) 

n N n,m  N n N n N 


(20) 


where  all  matrices  with  subscript  N indicate  network  quantities. 
Note  that  here  the  network  current  vectors  are  defined  by  referring 
to  the  tube  numbers.  An  alternate  way  is  to  refer  to  the  wavenumbers 
as  discussed  in  reference  (24) . 

Combining  Equations  (19)  and  (20) , we  obtain  one  form  of  the 
BLT  equation,  viz. 

„(3))„  - (P  „(s))  „)*(i^^^s))„  - (>^^®\s))„  (21) 

n,m  N n,m  N n N n N 

The  total  current  is  obtained  by  combining  Equations  (1)  and  (21) 


(I„ 

(s))  - 

((i  ) 

n 

N 

n,m 

where 

(1  J 

is  the 

n,m 

( 1 

1 

> > 

n,m 

0 

n,m  N 


n,m 


P ) )"^.(=/ 

n,m  If  n 


N 

(22) 


Kronecker  delta 
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SECTION  III 
EXAMPLES 


In  this  section,  an  example  of  a multiconductor  transmlsslon-llne 
network  Is  first  presented.  Numerical  results  for  a single  transmlsslon- 
llne  network  and  a more  general  multiconductor  network  are  then  given. 

1.  Example  for  a Mult Icon due tor  Transmission-Line  Network 

Consider  the  multiconductor  transmlsslon-llne  network  of 

Figure  la.  Let  us  label  the  tubes  and  junctions  as  In  the  graph  of 

Figure  lb.  The  junction  subgraphs  are  Illustrated  In  Figure  2.  In 

this  example,  we  first  demonstrate  how  to  obtain  the  matrices 

(P  (s))  (S„  „(s))„  and  (-3^®\s))^. 

n,m  N n,m  N n N 


(a)  The  propagation  matrix 

For  tube  T,  - the  current  propagation  matrix  (P  (s))_ 

1,2  2 ^1.2 

Is  given  by  Equation  (6)  with  (y  (s))„  given  by  Equation  (3). 

( ^3  ^ TQ  * 

The  tube  source  vector  ■’(s))!!’  Is  given  by  Equation  (7),  l.e., 

" ^1,2 


1,2 


2 ‘'n,m  ^1,2  " ^1,2 


ri  (’^n,m^.  -(expCy^  (s))  L 2)  • 
1,2  n,m  1,2 

^1,2  ‘'n,m  ^1,2  “ 


1.2 


where  )„  Is  formed  by  the  eigenvectors  of  the  matrix 

n,m  T^^2 

K Jr  '^K 

n,m  T^^2  ^1^2 

For  Tube  T«  - , we  have  (P  „(s))_  and  (y  (s))  given  by, 
2,3  n,m  ^23  ^ ^23 

respectively.  Equations  (6)  and  (3),  and  therA  Is  no  soufee  veetdr. 

Similar  procedures  are  applied  to  tubes  T^  ^ and  T2  ^ • 

The  network  propagation  matrix  (P^  network 

source  vector  can  be  readily  Identified  using  the  results  above,  and 

we  have 
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■ • - 


(b)  The  scattering  matrix 

For  junction  where  tube  ^ terminated  in  the  impedance , 

(z^  Equation  (18) 

n,m  '^1 

.W'j,  - «c  >1^  - '^c  >1,  > 

1 n,m  1 n,m  1,2  n,m  1 n,m  1,2 


For  junction  J2  , the  junction  connection  matrix  is 


n,m  2 


Tube  T, 


0 0 


Tube  T-  , Tube  T 

2,3  . 2,4 


-1  0 

0 -1 

0 0 


0 0 

0 -1 

-1  0 


The  corresponding  is 


n,m  2 


Tube  T, 


0 0 


Tube  T2^3  Tube  T^^^' 


-1  0 

0 -1 

0 1 

0 0 


0 0 

0 0 

0 -1 

-1  0 


The  matrix  of  all  characteristic  impedances  at  junction  Is  thus 


(q  (Z  (s)) 


n,m  c 


n,m  *^2 


-(Z  ) 

(0  ) 

n,m 


L(o  ) 

^ n,m 


(0  ) 

n,m 


(0  ) 

n,m 


(Z  ) (0  ) 

c ■^T^  > n,m 

n,m  ^2,3 


n,m 


(Z  ) 

C 

n,m  ^2,4 


The  scattering  matrix  is  thus  given  by  Equation  (16).  Hence,  collecting 
terms,  we  have 
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A '’'4, 3^ 

The  network  scattering  matrix  i®  obtained  by  suitably 

rearranging  the  above  junction  scattering  matrices  so  that  the  ordering 
of  the  components  of  the  Incident  and  reflected  currents  Is  the  same 
as  In  the  propagation  equation.  Thus,  this  equation  becomes 
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A 


IX 


Equations  (23)  and  (24)  contain  all  the  matrices  required  to 
solve  the  BLT  equation. 

2.  Numerical  Examples 

A general  computer  code,  QV7T,  has  been  developed  for  treating  this 
class  of  transmission-line  problems.  Both  time  harmonic  and  transient 
results  are  obtained,  the  latter  being  provided  by  a FFT  algorithm 
(ref.  25). 

To  demonstrate  that  the  above  approach  can  solve  transmission- 
line  problems  involving  closed  loops,  a simple  example  involving  only 
single  transmission  lines  is  first  presented.  The  use  of  single-line 
configuration  avoids  unnecessary  complications  Inherent  in  the  multi- 
conductor analysis  and  serves  to  Illustrate  typical  results.  Results 
for  a multiconductor  network  will  then  be  Illustrated. 

For  the  single-wire  example.  Figure  4 shows  a simple  five-j unction, 
five-tube  network,  which  contains  a single  loop  and  is  excited  by  a 
voltage  source.  The  characteristic  Impedances  of  each  line  are  75(2, 

g 

the  propagation  velocities  are  all  3x10  meters/sec  and  all  lengths 
are  1 meter.  All  Junctions  with  the  exception  of  junction  J5  have  zero 
Intrinsic  Impedances,  so  that  the  current  scattering  matrix  for  these 
junctions  depend  only  on  the  number  and  impedances  of  lines  connecting 
the  particular  junction.  Junction is  assumed  matched  with  a 75(2 
load.  A single  voltage  source  in  the  transmission  line  connecting 
junctions  and  J2  provides  the  excitation  for  the  network.  This  source 
is  located  midway  between  the  junctions  at  z ■ .5  meters,  and  is 
denoted  by  Vg  . 

Figure  5 shows  the  magnitude  of  the  voltage  at  the  75(2  load, 
denoted  by  , at  junction  , as  a function  of  frequency  for  a unit 
voltage  excitation.  The  various  peaks  observed  in  this  spectrum  are  due 
to  resonances  within  the  loop  structure,  as  well  as  resonances  on  the 
transmission  line  connecting  junctions  Jj  and  J^. 

25.  Cooley,  J.W.,  and  J.W.  Tukey,  "An  Algorithm  for  the  Machine 

Computation  of  Complex  Fourier  Series,"  Math.  Comput.,  Vol.  19, 
297-301,  April  1965. 


By  multiplying  this  impulse  response  spectrum  by  the  spectrum  of 
a step  function  voltage  and  performing  a numerical  Fourier  transformation, 
the  step  function  response  at  the  load  can  be  determined.  This  is 
illustrated  in  Figure  6.  Note  that  even  though  the  circulating  loop 
current  goes  to  infinity  as  time  progresses,  the  load  voltage  (and 
hence  load  current)  goes  to  zero,  as  required  by  a dc  analysis. 

A careful  examination  of  the  times  of  arrival  of  the  various  voltage 
pulses  provide  one  way  of  validating  the  correctness  of  the  analysis. 

Figure  7 shows  the  step  excited  load  voltage  at  junction  Jj 
for  the  same  network,  but  with  the  voltage  source  located  in  the  mid- 
point of  the  transmission  line  connecting  junctions  and  Jj . As  in  the 
previous  case,  the  times  of  arrival  between  the  various  waves  which 
comprise  the  total  transient  waveform  shown  in  this  figure  can  be 
used  to  verify  the  accuracy  of  the  calculation.  It  should  be  noted 
that  because  there  is  no  direct  current  connection  from  junction  J^,  J^, 

and  to  ground,  the  voltage  of  junction  J5  excited  by  a step  function 
excitation  Vg  , will  approach  zero  in  late  time,  due  to  the  lack  of 
current  flow  through  the  75^  load  at  the  junction  J^.  The  finite  slope  of 
the  vertical  lines  of  the  waveform,  most  evident  for  late  times,  is 
caused  by  the  fact  that  the  transient  waveform  is  sampled  at  discrete 
points.  In  the  ideal  case,  these  lines  should  have  an  infinite  slope. 

The  transient  response  of  the  load  voltage  due  to  the  two 
sources  previously  considered  exciting  the  transmission  line  at  the 
saem  time  is  illustrated  in  Figure  8.  Although  this  result  is  calculated 
directly  using  the  BLT  equation,  the  same  result  can  be  obtained  by 
superimposing  the  results  of  Figures  6 and  7. 

As  an  example  of  the  results  for  a more  complex  multiconductor 
network,  consider  the  transmission- line  configuration.  Illustrated  in 
Figure  0.  As  is  indicated  in  that  figure,  the  network  consists  of  four 
transmission-line  tubes.  Tube  T^^  2 ® three-wire  transmission  line  with 

tubes  T.,  , and  T-  , being  two -wire  lines . This  is  the  same  network  as  in 
Figure  1,  but  with  different  sources.  The  points  at  x/hich  the  network 
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Figure  7.  Step  excited  load  voltage  at  junction  Jj. 
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Figure  8.  Step  excited  voltage  at  junction  for  two 

driving  voltage  sources,  both  of  magnitude  Vg  . 
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response  is  calculated  are  also  f.ndicated.  The  linear  graph  for  this 

network  is  the  same  as  in  Figure  lb . 

For  this  example,  each  multiconductor  tube  length  is  assumed  to 

be  1 meter  long.  For  each  tube,  the  capacitive  coefficients  were  taken 

to  be  6,0x10  farad/meter  for  the  self  (diagonal)  terms  and 

-1.0x10  farad/meter  for  the  mutual  (off  diagonal)  terms.  The 

capacitive  coefficient  matrices  for  each  line  were  assumed  to  be 

highly  symmetric  such  that  these  two  numbers  completely  describe  the 

matrix.  All  propagation  velocities  on  the  multiconductor  lines  were 

0 

assumed  to  be  uniform  and  equal  to  3.0x10  m/s  so  that  the  inductance 

matrix  is  related  to  the  inverse  of  the  capacitive  coefficient  matrix. 

For  the  single-wire  transmission  line  in  tube  T , the  characteristic 

3,4 

line  impedance  was  taken  to  be  100  ohms.  That  impedance,  together  with 

0 

the  assumed  propagation  speed  of  3.0x10  m/s  determines  the  capacitance 
per  unit  length  of  that  line. 

The  terminating  impedance  network  at  junction  is  chosen  such 
that  the  line  is  perfectly  matched.  With  the  capacitive  coefficient 
chosen  as  above,  this  matched  load  is  described  by  an  impedance  matrix 
with  the  values  59.52  ohms  on  the  diagonal  and  11.90  ohms  off  the  diagonal. 
Other  terminations  for  the  transmission-line  tubes  are  indicated  in 
Figure  9.  Note  that  tube  T2  ^ has  its  open  circuit  at  junction  J^.  Similarly, 
the  second  wire  of  tube  T2  ^ open  circuited  at  junction  J^. 

In  this  example,  the  sources  exciting  the  transmission-line 
network  are  taken  to  be  current  sources.  A step-function  current 
source  of  magnitude  is  applied  at  time  t » 0 to  the  midpoint  of 

the  first  wire  in  tube  1.  Similarly,  there  are  two  current  sources  of 
magnitude  21^  applied  across  tube  T2  ^ ^ distance  of  .1  meter  from 

junction  2.  Finally,  tube  T2  ^ is  excited  by  a step  function  of  magnitude 
41^  occurring  at  time  t >•  0.  This  current  source  is  located  at  a 
distance  .4  meter  from  junction  12  and  excites  only  the  second  id.re  of 
tube  T2^4  . 

For  this  example,  the  voltage  response  at  three  load  points  is 
• considered.  These  load  points  are  Illustrated  in  Figure  9. 
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Figure  10  Illustrates  the  step  function,  current-excited  voltage 
across  load  point  1 as  a function  of  time.  This  quanitity  is  normalized 
with  respect  to  the  quantity  . Note  that  due  to  the  complexity 
of  the  transmission-line  network,  this  response  is  much  more  complicated 
than  responses  previously  seen. 

The  first  contribution  to  the  load  voltage  that  arrives  is 
due  to  a traveling  wave  from  the  source  on  tube  T^^  ^ load  point 

.5  meter  away.  This  corresponds  to  a time  of  1.66  nanoseconds.  The 
second  arrival  of  a wave  to  load  point  2 occurs  from  the  sources  on 
tubeT2  2 which  are  1.1  meters  away.  This  corresponds  to  a time  of 
3.66  nanoseconds. 

The  next  arrival  of  a wave  which  contributes  markedly  to  the 
voltage  response  at  load  point  1 is  at  9.66  nanoseconds.  This  contri- 
bution arises  from  a wave  traveling  along  the  first  wire  of  tube  T_  _ 

^ , J 

toward  junction  At  junction  there  is  a change  of  impedance  level 
of  the  line,  and  a fraction  of  the  wave  is  reflected,  eventually  returning 
to  load  point  1.  The  portion  of  the  wave  that  is  transmitted  through 
junction  continues  along  tube  T^  ^ and  suffers  complete  reflection 
at  load  4,  wehre  there  is  an  open-circuit  condition.  This  wave  eventually 
reflects  back  and  reaches  load  point  1 at  a time  of  t = 16.33  nanoseconds. 

All  of  these  times  may  be  verified  as  being  correct  from  the 
voltage  plots  presented  in  Figure  10.  Each  of  the  additional  break 
points  in  the  voltage  plots  of  Figure  10  can  be  attributed  to  arrival 
times  of  various  components  of  traveling  waves  on  this  transmission-line 
network. 

Figure  11  Illustrates  similar  results  for  the  voltage  at  load 
point  2.  Notice  that  the  first  time  of  arrival  is  at  t = .33  nanoseconds 
which  corresponds  to  the  distance  of  .1  meter. 

Figure  12  illustrates  the  step-function  current-excited  voltage 
response  at  load  3.  The  calculated  tum-on  time  for  this  case  is 
t “ 6.33  nanoseconds  and  agrees  well  with  the  results  of  Figure  12. 

An  additional  check  for  this  case  is  to  investigate  the  behavior  of 
the  current  at  load  point  3.  Since  this  point  is  open  circuit,  it 
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Figure  10.  Step  function,  current-excited  voltage  response  of 
network  In  Figure  9 at  load  point  1 . 


Figure  11.  Step  function,  current-excited  voltage  response 
of  network  In  Figure  9 at  load  point  2. 


is  required  that  the  current  be  zero  at  this  junction.  This  indeed  was 
found  to  be  the  case. 

Note  that  for  late  time,  the  normalized  voltages  at  load  points 
1,  2,  and  3 should  all  be  identical  since  they  are  on  the  same  wire 
with  no  intervening  Impedance  elements.  As  seen  from  Figures  10,  11 
and  12,  this  fact  is  verified  fay  the  computer  calculations. 
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SECTION  IV 
CONCLUSION 


This  report  has  described  and  Illustrated  a technique  for 
analyzing  lossless, multiconductor  transmlsslon-llne  networks.  This 
method  solves  for  currents  at  all  nodes  at  the  same  time,  as  opposed  to 
solving  the  problem  In  a junctlon-by-j unction  manner  as  In  the  conven- 
tional method.  The  present  approach  has  an  added  advantage  of  being 
able  to  treat  transmlsslon-llne  networks  that  contain  one  or  more  closed 
loops . 

The  ability  to  analyze  tansmlss Ion- line  networks  with  loops  Is 
Important,  especially  for  EMP  applications  where  there  Is  substantial 
energy  In  the  low  frequency  portion  of  the  spectrum.  Current  efforts 


Involve  extending  this  method  so  as  to  permit  the  analysis  of  general 
Interconnections  of  multiconductor  transmission  lines  which  are  excited 
not  by  lumped  voltage  and  current  soiirces,  but  by  distributed  sources 
which  are  provided  by  the  Incident  EMP  fields. 

The  analysis  of  a more  general  transmlsslon-llne  network  that 
may  Include  losses  will  be  reported  In  reference  (24) . This  upcoming 
report  also  discusses  In  greater  detail  the  general  theory  of  transmlsslon- 
llne  network  analysis. 
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